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(^ • Abstract 

The quantum statistics of the hght in the transparent medium with 
cubic nonhnearity is considered. Two types of transparent media are 
treated, namely, the cold transparent medium with a ground working level 
and the inversion-free medium with the lasing levels of the same popu- 
lation. The spectra of light fluctuation are found on the basis of both 
Scully-Lamb and Haken theories. The conditions for the use of effective 
Hamiltonian are determined. Basing on the exact solution of the Fokker- 
Plank equation for the Glauber-Sudarshan P-function the inversion-free 
medium with cubic nonhnearity is shown to be the source of spontaneous 
radiation with non-Gaussian statistics. 



Inroduction 



C^ ' Four-wave-mixing processes (FWM) are of particular interest for the study of 

nonclassical states of the light as predicted by Yuen and Shapiro Q and ob- 
served by Slusher et al. [||. The quantum theory of FWM is based on two 
general approaches. The first is the Effective Hamiltonian Approach (EHA) 
when a nonlinear medium is described by a c-number coupling constant in an 
effective Hamiltonian of interaction. The starting point of the second approach 
is the master equation for the field obtained by adiabatic elimination of the 
atomic variables. EHA has been used to consider the time evolution of the 
number of optical systems Q . Loudon et al have developed the quantum the- 
ory of propagation of the light and examined squeezing produced by self-phase 
modulation in fiber and in the nonlinear Saqnac interferometer with the help 
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of the effective Hamiltonian M. The effective Hamiltonians are usually in- 
troduced phenomenologically and occasionaly exact solution may be obtained. 
EHA seems to be attractive because it requires only susceptibilities to describe 
the nonlinear medium. For three-photon parametric phenomena the effective 
Hamiltonian can be obtained if the x medium is assumed to be transparent as 
it has been shown in . As to the case of x'^'^Hhe above mentioned conditions 
for the use of effective Hamiltonian construction would not be expected due to 
the existence of real transitions.lt is well known that the fluctuation-dissipation 
theorem (FDTh) does not exist in the case of four-index susceptibilities of a 
X^^^ medium Q. 

The second approach has been developed for two-level models: the Scully- 
Lamb model (SLM) for the excited-state atom Q and the Haken model (HM) 
for the ground-lower level atom ||§1 , |9) . Walls and Reid have shown that statistics 
of the light of nondegenerate FWM under linear theory over two weak modes 
is the same in both HM and SLM |^. For the single- mode FWM the master 
equations of HM and SLM are markedly different and the two models are in 
agreement only for the particular case of the linear theory B . 

The difference between HM and SLM arises not only from the type of the 
two-level atom. Formalism of the Haken theory is based on the multi-atom 
problem |10| . In contrast the starting points of the Scully-Lamb theory is the 
single-atom problem |jll|. In this view these theories are sometimes classified 
as the multi-atom and single-atom approaches p2[ . The problem is to find the 
most suitable approach [^, ||l3|. We believe that the HM and SLM describe 
an optical system in a quite different manner. It is an optical cavity with atoms 
inside for the HM while in the case of the SLM the atomic beam is injected into 
a cavity like in a maser or micromaser. It is clear that statistics of the light 
produced by the two systems may be substantially different. Here we use terms 
HM and LSM to take into account both the type of the approach (multi- or 
single-atom) and the type of the two-level atom. In the present paper we are 
aimed to apply our previous calculation within the framework of the Lax-Louisel 
model (LLM) to consider whether the effective Hamiltonians may be obtained 
for the x*-^-* medium and whether the statistics of the light is different for HM 
and LSM. 

We discuss in this paper an optical scheme including the cavity with non- 
linear medium inside and the classical injected signal. The x medium is 
described by the Lax-Louisell model which may be directly reduced to the HM 
and SLM. Two kinds of transparent medium are considered. In Section 1 spectra 
of fluctuation of the output light is calculated for the cold-transparent medium 
for which only the lower working level is populated and the frequency of the 
light is far from the atomic resonance. In Section 2 we consider the inversion- 
free medium when the population of both working levels are equal. The semi- 
classical Maxwell theory predicts no evolution of the light in the inversion-free 
medium. From this point of view the medium is transparent. Nevertheless it is 
a source of the light analogous with that of spontaneous radiation. With the use 
of the exact solution of the Fokker- Planck equation for the Glauber-Sudarshan 
P-function the statistics of the light is found to be non-Gaussian. 



1 The cold transparent medium 

1.1 Initial equations 

Let us consider a system consisting of N the similar three-level fixed atoms 
shown in Fig. la. It is the Lax-Louisell model when the incoherent pumping 
provides the populations /f and /| of working levels 1 and 2, respectively, 
under conditions of the field absence. The medium are placed into the high 
Q cavity (see Fig. lb) at w frequency, which can differ from the injected field 
frequency ul. The injected signal is considered to be monochromatic coherent 
wave. 




Fig. la. 
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Fig.lb. 

On the base of this model the master equation was derived by the adiabatic 
elimination of atomic variables in j^] where the field interaction with working 
transition was described by the Hamiltonian written in the dipole and rotating- 
field approximation. The resulted equation was the Fokker-Planck one for the 



Glauber-Sudarshan quasi-probability. 

dtP = (9a/ii + dl^f^2 + di^,fi3 + c.c.)P + R'pP, (1) 

where the coefficients fXi (correct to g'^) for the medium with cubic nonlinearity 
are: 



A^i = ^^«(i-/3i"n(/r-/i) 

+ ^(/|-/f)) + f^/|] 

^3 = -^^^{^(/2 - /r){^/2 + ^(1 - /r) (2) 

+/I - /i] + /3|aP^T^(/l - /i') + /3|apf /I - /|} 

i 1 + i 2 ^7 

(5 = 7 + iA 

A = UJ21 — UJ 

7 and ^21 are the constant of the transverse relaxation and the frequency of 
the atomic transition respectively, g — dTi~'^ {fiuj /2V eoY^"^ is the couphng con- 
stant, d is the dipole moment of the working transition, and (3 — 2g'^^{Vi + 
r2)(rir2|(5p)~"'^ is the saturation parameter. The relaxation constants Fiand 
r2 are formed from the set of initial constants: 

r2 = (7i + 72) + (7T + 7i)(Ai + A2 + 72)(Ai + A2 + 7i)~' 

7;=7fe(l-Afe(Ai+A2)-i) 
fc= 1,2. 

The constant Fi can be obtained from r2 by replacing . 7| ^ 7x, 7x ^ 7|, 
2 -> 1, 1 -> 2. The term 

Rp = da{-Ca — flg exp(i(a; — ujl)^)) + c.c. 

simulates the input and output of radiation from the cavity [pl,||I3. Here 
C = an + Cout, Cin(out) = cT'"^°"*^ L-\ c is the light velocity, T™(°"*) is 
the transmission coefficient of the input and output mirrors, L is the cavity 
perimeter, and aa is the amplitude of the injected field in inside the cavity. This 
amplitude is considered to be real. 

Equation (1) describes the linear and nonlinear processes of interaction. One 
can introduce the macroscopic susceptibilities for these processes. 

X^^H^) = Ng^{LJ2i-uj-i-f)-^ 

X^^Hu;) = Px^'Hu;) (3) 



The Maxwell equation for the complex amplitude of the field inside the cavity 
z = J (PaaP{a,t), follows from Eq.(l) Accounting for Eq.(3) this equation can 
be written as 

C 

dtz = {{iKi - K2)z + {-ixi + X2)z\zf){f^ ~ f^) - — z + aa eyip{i{uj - ujL)t) (4) 

where the real and imaginary parts of the susceptibilities are separated out 

X = Kl + iK2 

X^^^ = Xi+ *X2 

Since Eq. (1) is written correct up to g"^, the dimensionless intensity of the field 
cannot be large. 

f3\z\', P{\a\^) «1 (5) 

We consider the cold transparent medium assuming detuning x = A/7to be 
large. 

|x| » 1 (6) 

Under these conditions the real and imaginary parts of the susceptibilities are 
of essentially different order. 

K2 = Ki/x, X\=Xll^ 

Therefore the linear and nonlinear dispersion described by the susceptibilities 
Aviand xi have a dominant role in a transparent medium where /| = 0. The 
processes of absorption can be neglected setting ki,X2 = 0. However the last- 
mentioned fact means that the main losses of the field in the cavity are connected 
only to the radiation output, i.e., ki-, Xi\z^ ^ C/2. 

As a result for the cold (/| = 0) transparent medium we have: 

A*i = [-iKia + iXia|ap]/i 

M2 = aV(i + ri/r2)-Vr(-f(i-/r) + 

+*(/f - 1 - /rri/27)) (7) 

M3 = |apxi(i + ri/r2)-Vi^(i-/o-r 



2 

The Lax-Louisell model allows a direct transition to HM and SLM. In Q the 
necessary conditions were formulated which can be reduced to the set of in- 
equalities between constants 71,7^ and 71,72- 
So, putting in Eq.(7) 

Ti = r2=7T+7i=7||, /r + /l = l 
7 = 7-L 
711 /27± = /c 



we obtain the Fokker-Planck equations for HM. 

dtP = (a„(-iKia + ixia\a\^) + ^^(-^Xia'/c) + c.c.)P + R'pP, (8) 

In this equation the various types of relaxation are described < /c < 1, the 
fc characterizes the difference between longitudinal and transversal relaxation 
velocities. The case when fc = l describes the radiation decay only. We note 
that Eq.(8) agrees with equation obtained in WM. 
Putting in Eq.(7) 

Tfe = 7fc 
ifkf = 0, /I^O, k = l,2 (9) 

we obtain the SLM equation. 

dtP = f^{dai-iKa + ixa\a\'^) 

+^L(-Xi(l + 7i/72)-'«'(* - x/2)) (10) 

+dL' (Xi(l + li/l2)-'\a\^x/2) + c.c.)P + R'pP 

The conditions (9) include the requirement for weak pumping (/^)^ = 0,/^ 7^ 0, 
when the terms proportional to (/^)^ are absent. Formally, it is the consequence 
of the single-atom approach of the SLM. Physically, SLM describes optical sys- 
tems in which the atoms are injected into the cavity being in a prepared state 
and rapidly escaped from the interaction zone. The atoms output from the cav- 
ity is modeled by their output from the working levels. On one hand, it requires 
an introduction of two-level systems with excited levels, and, on the other hand, 
it allows the use of the single-atom approach. 

The Lax-Louisell model discussed here a priori describes the atom with two 
excited working levels. As it follows from Eq.(7) at 

fl = i, rir2((ri + r2)7)-' = /c (n) 

the Fokker-Planck equations in the Haken and the Lax-Louisell models are the 
same. It means that the statistics of the field formed by two-level atoms with 
ground state and that of the field formed by atoms with two excited levels is 
the same. Therefore under conditions (11) the differences between LSM and 
LLM mainly depend upon the type of physical system rather than the type of 
two-level atoms. We note that these differences disappear when the linear over 
interaction theory is used. This theory is used as an example in the case of 
the linear medium or in the case of multi-photon interaction described by the 
effective Hamiltonian of the field in the linear approximation |l4| . 

1.2 The effective Hamiltonian 

The effective Hamiltonian introduced phcnomcnologically is often used to de- 
scribe the four-wave mixing in the transparent medium with cubic nonlinearity. 



The Hamiltonian is 

Heff = hka+\^, (12) 

where k is the couphng constant. The exact equation for the Glauber-Sudarshan 
quasi-probabihty written on the base of Eq.(12) has the Fokker-Plank equation 
form 

dtP = (a„(ito|an + dl^{~'-ka^) + c.c.)P. 

For the cold transparent medium described by the Lax-Louisell model the ef- 
fective Hamiltonian appears solely in the case when 

/r = 1, 

rir2((ri + r2)7)"' = i. (i3) 

In the HM it means that /c = 1, i.e., the decay is of the purely radiation 
character. Under conditions of Eq.(13) the equation for the electromagnetic 
field density matrix p following from Eq.(7) or Eq.(8) takes the form 

dtp = -irr^[Hi +H2,p] + RfP 

H2 = \nxia+\'' (14) 

where the operator Rpp in P-representation corresponds to the R'pp operator. 
In contrast to the effective Hamiltonian approach (EHA) based on Eq.(12) the 
effect of the linear dispersion described by the Hamiltonian Hi is taken into 
account in Eq.l4. 

Physically speaking the linear dispersion tends to the change of mode fre- 
quency oj ^ uj — Ki. It results the additional detuning between the injected 
field frequency and the cavity mode. The role of this factor was discussed in 
[^ and, as will be shown, it substantially affects the fluctuation dynamics and 
the noises. If the injected signal frequency is such that detuning is absent, i.e., 
ojl — (w — ki) = 0, the linear dispersion is not important in the transparent 
medium. In such a case the description by EHA with k = Rex^^' coincides with 
that by Eq.(14). The effective Hamiltonian principally does not arise in the 
SLM due of the terms proportional to x. 

1.3 The Langevin equations 

The condition of the small fluctuations allowing the linearization of the Fokker- 
Planck equation (FPE) is used as a basic approximation for the calculation of 
observed values. Basing on the linearized equation we write the Langevin equa- 
tions which are convenient to obtained the fluctuation spectrum. We introduce 
polar coordinates 

aa* = I, ip = (l/2i) ln(a/Q;*) 



The fluctuations of intensity and phase around its semi-classical values are as- 
sumed to be small 



I=U + S £< U 

(fi = ipo+'^ ip '^fo 



(15) 



where U = |zp and (po = a,rgz are the steady-state solutions of the semi-classical 
equations (4), which for transparent medium take the form : 



dtU =-CU + 2anVUcosLpo = 
dt(p = Ki - xiU + aosmipa/VU = 



(16) 



In Eqs.( 16 ) we assumed the phase of the injected field to be equal to zero and 
ujl = oj- As a result we obtain the Langevin equation 



a 



^ 



-A, Ai2 
^21, -A 



l]+D^'Mt) 



(17) 



where ri(t) is the two-dimensional Winer process and the elements of the drift 
A and diffusion D matrix for the three considered models are listed in Table 1. 
Note that the linear dispersion affects only the fluctuations dynamics, and in 
SLM the source of the phase noise appears. 



Table 1. 



The phase and the intensity fluctuation are coupled dynamically A12, A21 ^ 
and statistically D^^p ^ 0. This fact leads to two circumstances. First, both 
the intensity and phase fluctuations is produced by two sources Demand D^^. 
Second, under steady-state the fluctuations being attenuated oscillate with the 
frequency 12. To obtain 17 we write closed equations for the averaged values. 
These equations follow from Eq.(17). 





EHA 


HM 


SLM 


A 


C/2 


A12 


2kU'^ 


2xiU-' - 2Uki 


A21 




-iXi + i«i/C/ 


0,4, 


-kU 


-XiUfc 


-XiU 


Dtp^p 





l^Xi 



{d^t + 2Adt + A^ + ^'^)(j 



n 





-A12A21 



(18) 



In HM and SLM 



Vi' 



k2(1 -/?[/) (1-3/3C/), 



where the dimensionless intensity (311 <C 1 since the initial equation (1) is correct 
up to (7^. Let 



PU < 1/3 



(19) 



Then Q^ = K^and fi takes a new meaning of the damped harmonic oscillator 
frequency described by Eq.(18). As a result under steady-state the fluctuations 
of phase and intensity damp with velocity C/2 and oscillate at the frequency fl. 
In HM and SLM under conditions (19) the dynamics of the fluctuations are 
defined solely by the linear dispersion. 

Ai2 = -2Uki, A21 = -Ki/U 

where the expression for U follows from Eqs.(16) under conditions (19). 

U = al{A^ + Kl)-' (20) 

The output power P = hujcUT°'^^L^^ will be connected to the power of the 
injected light Pq with the use of Eq.(20) as 

1 4e 

P = P^- 



1 + ^2 (l + g)2 

where 

t = 2ki/C, e = r"/T°"* (21) 

1.4 The fluctuations spectra 

For the homodync balanced detection, when the signal is mixed with the field of 
local oscillator, the normal-ordered correlation functions of quadrature operator 

X(t,e) ==a+(t)exp(ie) + /i.c. (22) 

are measured. The spectrum of the photocurrent i^ \lij) or the noise power 
spectrum can be presented in the form: 

g{LU,e) = C™* / dTe'^^{:X{t,e)X{t + T,e):) (23) 



where 77 is the quantum efficiency of the recording scheme, the points denote 
the normal ordering, and scale is determined so that the unit corresponds to 
the shot noise level. The value 5(0;, 0)can be calculated by introduction of the 
Langevin variable 

X{t,Q) = a*{t)exp{ie) + c.c. (24) 

which is c- number analog of the operator (22). So, it is easy to show that 

/+00 
dTe"^''{X{t,e)X{t + T,e)) (25) 

'CO 



The Eq. (25) is true when the Langevin equations are written from the FPE 
for the Glauber P-function or, as an example, the generalized P-representation. 
With the use of the small fluctuations (15) 

X{t, e) = 2VUcos(po + 2VUsm{e - ipo)'^(t) + -^ cos(e - (po)£{t) 



Therefore the case when the phase of the local oscillator field coincides with 
that of the signal, i.e. O = i^Oj corresponds to the measurement of spectrum 
of the intensity fluctuations, while when O — (po = § the phase fluctuations are 
measured. Solving the Langevin equation (17) we obtain 



ff(c.,e) 



= 2C°"*[('-^ 



■{D^^[{^lA + vAi2f + M^tj2] 
+2D^^[{nA + uAi2){plA2i + vA) 



jiiyu) 



where ^ = y/WJ sin (6 — (po), v = —= cos (9 — (^o), and 5*1,2 are the roots of 
the equation {S + Af' — A12A21 = 0. Introducing the scaled frequency ui = uj/A, 
t = Ki/A we represent the obtained spectrum in the form 



g{Lo,e) = 



4 W{Q-ipo)-t^ + V{e-ipa)-u;^ 



1 + e 



((1-2 + 1 - i2)2 + 4^2 





EHA 


HM 


LSM 


W{0) 


-2/3 


WUfc 


2/3C/(l + f ) 


V{0) 





wi^) 


2 


-2/3C//e 


_2/3t/(l-^) 


V{^) 





PUxt 


t 


VSkU/A 


«iM 



For the measurement of the intensity and the phase fluctuations the coefficients 
are presented in Table 2. For SLM the case 71 = 72, /i = lis taken as before. 



Table 2. 



When EHM and HM are compared one can sec that the effect of dispersion is 
diametrically opposite behavior both for the phase and amplitude fluctuations. 
However on the basis of both approaches compaction arises leading to a decrease 
in the shot noise by a factor of not more than 3. Indeed the noise reduction 
may be obtained in nonzero frequency range. This behavior is understood by 
suggestion that the dispersion rotates the squeezing ellipse. 

To determine the degree of squeezing for the optical scheme under consid- 
eration we need to choose the optimal phase of the local oscillator Q — Qq. 
To this end we represent g(u;, O) in the form g = A + B sin 26 + C cos 20 and 
determine the 69 value from the condition tan26o = B/C. As a result the 
fluctuation spectrum takes a form 



5(^,60) 



1 + e 



[{uj^ + l-ty+At^]-'G{uj) 
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where 



GEHAiuj) - |(2t ± (4t^ + 3(1 -t^+ Co^ff'^) 

GHM{uo)^±2f,pu\t\-M 

GLSAiioj) - 2l3Uxt[x{l + cj2 + f) ± M] 

M = [(1 + t2)-l[(l +Q^+ ff + f{l +f- (I-2)2]]l/2 



In the zero-frequency range 



4 t 



gEHA{^,Q)^^j^-^^{2t±./WTW^W) 






9hm{0, e) = ±4/e/3£/ , ' ', ^ sguee^ing - 2/3 



4 t 

5lsm(0, O) = ±—l3U— weaksqueezing 

X t-^ + 1 

For EHM the noise can be suppressed not only in the zero-frequency range. 
Choosing the local oscillator field phase using the condition tan0 = (-s/Si)"^, 
wc find 



gEHAi^) 



8Q^ i 



2 



At the frequency ujg = 1 + t'^ we find the minimum 



grain \^e) 

The maximum drop is equal 2/3. 



2t 



2 



C(l + 3t2) 



2 Inversion- free transparent media 

2.1 Initial equations 

Let the pumping produce equal populations of working levels fl — fi — fo- As 
it seen from FPE (1), in such a medium there is no polarization proportional to 
inversion. In this sense the inversion-free medium is transparent. Nevertheless 
the field can evaluate in such a media, especially in nonlinear one. At the same 
time it is seen directly from Eq.(l) that the contribution of second derivatives 
formed from the terms of /i2 and /is remains depending solely on the population 
of upper level /| so the field can nevertheless evolute in such mediua, espe- 
cially in nonlinear ones. Therefore the field in a inversion- free medium has a 
spontaneous radiation nature. 
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However the Eq.(l) does not permit the correct consideration of the field in a 
inversion-free medium because of the procedure of its derivation. The reason is 
that the noncommutativity of differential operators is neglected when the second 
derivatives are obtained. It has been shown in M that the noncommutativity 
contributes to the first derivatives in FPE that correspond to the semiclassical 
Maxwell's equations. If ff ^ /| this contributions can be neglected. However, 
for inversion-free medium the first derivatives are produced by the noncommu- 
tativity only that why it must be taken into account. Using the formalism 
from [^ it is possible to obtain the master equation for the field with due re- 
gard for noncommutativity of differential operators. The Glauber-Sudarshan 
quasi-probability equation correct to g^ for LLM has the form 

I ^ ro 1 II "i iM |2 1 , 2/o(r2 - Ti) - r2,, , , 
-[2K2 + X2((a; -l)|a| -1 + T^ , T^ )] (26) 



dada*' -^"^ '' ' Ti+Ts 

^' 9 h . ^2 - Ti , . (/o-i)(r2-ri)^ 



da*da^ -^ ' ' T1+T2 Ti 

^' r2x2/0Tr^ +c.c)P 



da'^da*'^ -^ •"Ti+r2 

The polarization proportional to cubic susceptibility X2 appearing in inversion- 
free medium leads to evolution of the complex field amplitude (a) = J (PaaP{a, t) 
So, according to Eq.(26) 

5t(a)=-^/o[X2(3 + x2)-2zxi](a) 

The field amplitude damping and its steady-state value equals zero because 
X2 = 0. From the physicall point of view this result is related to the diffusion 
of the spontaneous radiation phase. Equation (26) can be approximated by the 
following FPE 

dtP = \h{^ax2[i + x^-2ix\-^a'n2{l + x^) (27) 

2 oa oa'^ 

^' -[2« + X2{{x^ ~ l)|ap - I)] + c.c)P + R'pP 



dada*' ^ " " ' 2' 

where we set Fi = r2 and add the term RpP assuming that the medium is 
placed in the cavity with the injected signal on. The FPE obtained can be 
solved exactly. 

2.2 The exact solutions for free-inversion transparent medium 

Let us represent the FPE (27) in the form 

dtP^~da{ka + ao)+daa'^a^Aaa (28) 

+aL.(6+|a|'A„a.) + c.c)P 
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where Uq and b are the real constants. To obtain the solutions the corresponding 
Langevin equations are used 

dta — ka + uq + fa 

dta* ^ (ka + ao)* + U' (29) 

where nonzero correlation functions for random forces have the form 

{U{t)fait + T))=2{a^)AaaS{T) 

(/a(t)/a* {t + r)) = 2(6 + (|a|2)A„«.(5(r) 

The formulated Langevin approach based on the FPE for the Glauber-Sudarshan 
quasi- probability was clasified in ||lj| as a "naive" one, because, for instance, 
the variables aand a* cannot be interpreted as independent ones. Here we do 
not discuss neither the interpretation problems nor the mathematical aspects of 
stochastic differential equation theory, but only notice the following. First, the 
Langevin approach is in full accordance with the FPE (27) as far as concerned 
the calculation of observed values and based on Itto interpretation of stochas- 
tic integral. Second, all the results obtained below can be readily reproduced 
using, for example, the generalized quasi-probability P{a,a^), where aand a^ 
are independent variables owing to the simple correspondence 

P{aTa*) < — >P(a, a^) a,a* < — > a^a^ 

Solving Eq.(29) we find 

ait) = zit) + j dse''^'-"'^ fa{s) 
Jo 

z{t) ^ a{0)e''' + ^{e"' - 1) (30) 

The expression for steady-state of complex amplitude of the field in the cavity 
can be consequently obtained. 

(a) == -flo/fc 
i?efc<0 

Using Eq.(30) we find the correlation function 

Y{T,Q) = {x{t,e)x{t + T,e)), T>o 

X{t, e) = a* exp(ze) + c.c. 

To do this we need in average values of the type {a{t)a{t+T)), {a* {t)a{t+T))a,nd 
so on. For example, using Eq.(30) we obtain: 

{a{t)a{t + T)) =^ z{t)z{t + T) 

+e''*[{a^{t))-z\t) + e^''' [ dsi [ ds2 

Jo Jo 
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Due to the delta-function properties this integral is equal to zero, consequently, 
at i ^ oo 

{a{t)a{t + r)) = (a)^ + e''-{{a^) - (a)^) 

As a result at t ^ oo the expression for Y takes the form 

Y{t,Q) ^ {a)(a*) + e-^'^{af 
+e'^'{{\a\')-{a){a*) 
+e^''^i{a^)-{af))+c.c. (31) 

The expression for the Fourier transformation is written as 

/•CO 

Y{uj, 6) = 2Re / dTe''^''Y{T, 9) 
Jo 

Y{lu, Q) = ^[(CJ^ + 1 - T^f + 4T2]-i . 
■{{S + Ren){l +il)^ + T^) + TImfi{Oj^ - 1 - T^)) (32) 

where 

T ^ Im k/A 

UJ = LO/A 

A = -Rek 

M = exp(-2ie)((a2)- (a)^) 

It should be noted that Eq.(31) is valid for any diffusion coefhcient of FPE (28). 
Therefore the main problem reduces to the determination of the steady-state 
variance S and /x. 

Using Eq.(28) or the Langevin equations (29) the closed equations for S and 
II can be obtained, which can be solved exactly. As a result, at 

i?eA: < 

Re{k + Aaa') < 
Re{k + Aac) < 

we obtain the necessary steady-state solutions: 

(l«l') = {-^{k + k*)-2b){k + k*+2A^^.)-' 
|fc| 

^ ^ -2(^A„„. +6)(fc + r+2A„„.)"' 



-e 



i2 A 

2ie "0 ^^aa 



|fc|2 k + Ao 



(33) 



It is interesting to note the following fact. The Fokker-Planck equations (8) 
and (9) for the SLM and the HM models reduce to Eq.(28) if the nonlinear 
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dispersion is neglected. For example, assuming that k ^ Xi|q^Pi i-^-, /3|ap ^ 1 
we obtain in Eq.(28) 

k = -C/2 + zKi 

Aaa - (-^Xlfc + Xlx/4)fi 

Aaa* = /ra;xi/4 
6 = 

where fl ^ 1, x — should be set for the Haken model and a; ^ 0, /c == 1 for the 
Scully-Lamb one. The saturation parameter /3 can be considered in transparent 
medium as a small one (in the sense of the conditions (6)). This permits the 
neglect of Aa^and Aaa* values in denominators of Eq.(33) when calculating S 
and fi. This approximation properly implies a linearization of the FPE near 
the steady-state (a)obtained according to Eq.(30). As a result the solutions 
obtained in Section 1 both for the Haken and the Scully-Lamb models can be 
found using the presented formulas. 

2.3 Fluctuation spectra 

The term |/o^ax2(— 2«a;)in FPE (27) corresponds to the dispersion leading 
to the change of cavity mode frequency if w ^ cj' — /0X2 ■ 2;. However, if the 
injected signal frequency is assumed to be ujl = ^' the dispersion plays no role. 
For this case the Eqs.(33) take the form 

(l«P) = ('7o + ^(l + i/3?o(3 + x2))i? 






S = (90 - :^go(l - a;2))i? (34) 



M = -|^go(l+2;2)i?exp(-2je) 
R = (1 + 2/3(70 



,-1 



where the assumption go = 2/0K2/C' is used and 1 — 3/3/4 w IBecause the value 
k = —^[1 -I- ^/3(7o(3 -I- x^)] is real the expression for spectrum Y{ij^ Q) obtained 
according to Eq. (32) takes the simple form 

y(.,e)^ 8(^ + ^eM) 



C[l + i/3go(3 + x2)](l+<:s2) 
Both the maximum and minimum values of Y is obtained at = -^ (max) 

S + Ren = go(l + x^l3al/k'^){l + 2/3go)-i 
and = (min) 
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S + Refi = qo{l - I3al/e)il + 2/3go)"' 

. For the case = 0, which corresponds to the measurement of amphtude 
fluctuations, the value S+Reii > 0. It follows from the fact that initial equation 
is correct to the order oi g*i.e., only the cubic nonlinearity is taken into account. 
Therefore it is necessary to consider the field as not very strong one, i.e., 

pal/k^ ^f3{a)(a*) <1, /S < \a\^) < 1 

As a result the noise is found to be greater that the shot one. Physically it is 
clear because the inversion-free medium is the source of spontaneous radiation, 
the statistics of which, however, can be non-Gaussian. 

2.4 Generation of light by the inversion- free medium 

Let the injected signal and detuning be equal to zero. In such a medium the 
field with average number of photons 

(n) = (7o(l + 2/3go)"' 
qo - 2foK2/C (35) 

and the Lorentzian profile of spectral line with halfwidth 

is generated. Notice that the threshold of the generation is absent, and dimen- 
sionless intensity is limited: 

/3(n) < 1/2 

Consider the statistics of the generated light. Using Eq.(27) the Mandel param- 
eter can be obtained: 

For the linear inversion- free medium /3 = 0, ^ = (n)i.e., the field statistics is of 
the Gaussian type. By contrast in the case under consideration |(n) < C < {^■ 
It means that the function of the photons number distribution is narrower than 
the Gaussian one. Therefore the spontaneous radiation statistics of the non- 
linear inversion-free medium is of essentially non-Gaussian nature. 



3 Conclusion 

The presented study of quantum statistics of single mode light in transparent 
x(3) medium allows us to compare the intensity and phase fluctuation spectra 
obtained on the basis of both Haken and Scully-Lamb models. The performed 
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calculations clearly show how the light statistics in transparent x(3) medium 
depends on the model medium excitation. It is seen that the light in Scully- 
Lamb model is more noisy than in the Haken one. Physically, it is rather 
evidently since in the SLM the medium subsystem is open. The atoms are 
injected into the cavity and leave the interaction zone rapidly. As a result, the 
SLM the medium excitation has additional source of noise whith relation to the 
Haken model. It is shown in this paper that there are no any cases when EHA is 
valid in framework of Scully-Lamb model. But for the cold transparent medium 
Haken model with the lower working level to be ground state the traditional 
effective Hamiltonian of interaction was derived in particular case of absence of 
the collision-induced phase decay of atoms. 

In the interesting case of the inversion-free x(3) medium presented study 
leads to description of quantum evolution of light based on reduction of the 
Fokker-Plank equation from Q to equation (20) presented paper. The obtained 
FPA has an exact solution. The statistics of generated light in this case is 
nongaussian that is photon number distribution is narrower than in the Gaussian 
statistics. At the same time there is no squeezing. Is this state of generated 
light quantum because this state and its evolution appear only in the framework 
of quantum description? It necessary that futhcr study of the problem should 
be performed. 
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